Time-dependent evolution of two coupled Luttinger liquids 
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We consider two disconnected Luttinger liquids which are coupled at t — through chiral density- 
density interactions. Both for t < and t > the system is exactly solvable by means of bosonization 
and this allows to evaluate analytically the time-dependence of correlation functions. We find that 
in the long-time limit the critical exponent governing the one-particle correlation function differs 
from the exponent dictated by the equilibrium ground state of the coupled system. We also discuss 
how this reflects on some physical quantities which are accessible to real experiments. 



Time-dependent quantum systems are attracting much 
the attention during the last years. This is because the 
dynamical response of nanoscopic devices is nowadays 
accessible in many experiments with high relevance in 
practical applications such as quantum computing and 
single-electron transport 0. On one side we mention the 
wide investigation of out-of-equilibrium phenomena in 
transport experiments, where time-dependent transient 
currents are measured at picosecond time-scales On 
the other side, recent experiments in ultracold atoms 
confined in optical lattices have shown that it is possi- 
ble to time-tune the strength of the interactions in both 
bosonic .3j and fermionic|J] systems, by using the so called 
Feshbach resonance. In these cases it arises the intriguing 
question of what is the steady state to which the initial 
ground state relaxes top\ 6]. 

In the light of this challenging physics, the theoretical 
investigation of the time-dependent evolution of many- 
body interacting systems deserves special attention. Re- 
cently, an exact formulation of time-dependent transport 
in electronic systems 7] was derived in the framework 
of time-dependent-density-functional-theory 0. Unfor- 
tunately its implementation to strongly interacting sys- 
tems results hard. An efficient numerical tool to address 
the problem of the electron-electron (e-e) interactions 
is the time-dependent-density-matrix-renormalization- 
group00]0|5 which turns out to have excellent perfor- 
mance in one-dimensional (ID) systems. For what con- 
cerns electronic systems with time-dependent parameters 
but not involving charge transport, in a recent paper 
Cazalilla carried out the exact evolution of a noninteract- 
ing ID system following a sudden switch-on of Luttinger 
liquid e-e interaction[12j. Remarkably, it is found that 
the initial ground state reaches a stationary state only if 
the system is infinite-sized and that its asymptotic be- 
havior differs from the one of the interacting equilibrium 
systenvfjf . 

In the present paper we extend this analysis to two 
coupled interacting ID systems. Our system consists of 
two isolated Luttinger liquids which are connected at t = 



through interactions between chiral electron densities 
located at different liquids. 

The Luttinger liquid is the prototype of interacting 
electrons in ID and it is governed by the so called 
Tomonaga-Luttingcr Hamiltonian. In this model, the 
electrons have linear dispersion relation around positive 
(Right) and negative (Left) Fermi points located at ±kp 
and the e-e interactions act only between Right /Left elec- 
tron densities. This means that only the component with 
momentum transfer p ~ of the Coulomb repulsion is re- 
tained, while the p ~ 2k}? component is assumed to be 
negligible. The model is exactly solvable by means of 
bosonization technique, which allows to write the elec- 
tron Hamiltonian in terms of boson operators &'s. 

In the following we shall adopt a similar notation as in 
Ref. [l3 

For t < the system consists of two disconnected 
identical interacting electron liquids described by the 
bosonized Hamiltonian 0] Hq = Hi + H2 with 

H * = \Y. ^F+9 ( '\q))\q\{b\{q)bM) + h{q)b\{q)) 

-9t\q){b\{q)b\{-q)+b i {q)b i {-q))] (1) 

where i = 1,2, \bi{q),b-{q')} — 5i^5 q , q ', vp is the Fermi 

velocity, \q) is the interaction parameter between 
Right-Right (positive q) and Left-Left (negative q) elec- 
tron densities, while g\ (q) is the interaction parameter 
between Left-Right densities. 

The system is diagonalized by the well-known Bogo- 
lioubov transformation: 

bi(q) = cosh (p(q)bi(q) + sinh <p(q)b\(-q) , 

b\{q) = smh^(q)b l (-q)+coship(q)bl(q) (2) 

with 

tanh2^) = gf\q)/[v F + <?f (<?)] , (3) 
and renormalized velocity 

v(q) = ■\]{vF+9t\q)?-{9 < iXq)Y- (4) 
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For t > the chiral density-density interactions be- 
tween the two Luttinger liquids are switched on and the 
system is governed by the total Hamiltonian 



H to t = H + 0(i)H 12 



(5) 



where 



H12 = 



bUq)h(q) 



bl(q)bl(q)+b 2 (q)b\(q) 

^ ) (9)(&l(?)4(-9) + M<zM 

bUq)b\(-q) + b2(q)b 1 (-q))], 



-q) 



(6) 



where g^ (q) ( g^ (<z) ) is the interaction parameter be- 
tween electron densities of the same (opposite) chirality 
in different liquids. Physically H\ 2 could represent the 
long-range component of the Coulomb repulsion felt by 
electrons located in two (quasi) ID metallic systems close 
to each other. Here we are considering a general situa- 
tion where g£ (q) ^ g^\q) but in real systems like car- 
bon nanotubes, they coincide (see below). It is worth to 
remark that in the present model the coupling between 
the two liquids does not involve any inter-liquid electron 
tunneling. 

H to t is again diagonalized by introducing symmetrized 
and antisymmetrized operators |l4): 



s(q) = [h(q)+b 2 (q)]/V2, 
a(q) = [b 1 (q)-b 2 {q)]/V2. 



(7) 



In terms of them H to t decouples in two independent non- 
equivalent Luttinger liquids: 



H s + H a 



= ~ T,^ + 9t ) (q) + 9 i 2 ) )\q\^(qHq) + -<q>Hq)) 



q^O 

,(2)/ s , J2), 



H9r(q)+92 2 \q))^(q>H-q) + S (q)s(-q))} 
4 E^ + - 9 i 2 ) )\q\{a\q)a{q) + a{q)a){q)) 



J 2 ), 



-(9\ 2 \q) g\ l \q))(a\ q )a\- q ) + a(q)a(-q))] . (8) 
Finally, the diagonal form of H to t is obtained in terms of 



Hq) 

a(q) 
aUq) 



cosh ip s (q)s(q) + shahip s (q)s^ (-q) , 
sinhip s (q)s(-q) + cosh ip s (q) s f (q) , 
coship a (q)a(q) + sinh<£ a (q)a t (-g) , 
sinh ip a (q)a(-q) + cosh ip a (q)a\q) , 



(9) 



where 



tanh2p„ ( g ) = [g{ 2 \q)±g! 2 2 \q)}/[v F +gi 4 \q)±g! 2 i \q)} , 

(10) 



and renormalized velocities 



VsAq) = \l (v F + sf (q) ± g ( 2 A) {q)Y- ( 9 i 2) (q) ± 5 f (q)) 2 

(11) 

Now we are in the position to evaluate the equal-time 
one-particle correlation function for t > 0, defined as 



Gf OM) = {e iH ^ in (x)^(0)e 



-Hit. 



H 



(12) 



where i = 1 , 2 labels the two Luttinger liquids decoupled 
at t < 0, the subscript 7 = R, L indicates the Right/Left 
character of the non- interacting fermion fields ip and tp^ , 
and (. . .) h is the zero-temperature average in the en- 
semble Ho. Without any loss of generality we shall focus 
on G$(jc,t). 

The bosonization technique allows to calculate Gr 
with logarithmic accurancy, which is exact for dis- 
tances/times much longer than the typical range of the 
interactions. The key point of bosonization is that it 
is possible to express the fermion field in terms of boson 
fields. For instance the Right-mover fermion field is given 
by the following expression: 



(13) 



where a is a short-distance cutoff, proportional to the 
lattice spacing, t]r is an anticommuting Klein factor and 

q>0 ^ ' 



ipo,R + 2-kxN r /L , 



(14) 



where Nr is the total number of Right-electrons, 
[<fio,R, Nr] = i and L is the length of the system. 

Thus the computation of the correlation function re- 
duces to: 
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i*l,ii(a:,t)p-i*l,H(0,t)i 



Ho 



JL e -H<(* i .«( a! ' t )-* i . ii ( 3! ' t )) 2 > ir D-[*i.«( a! ' t ).*i.«( o . :t )]}(i5) 

where $i tR (x,t) = e iHtott $ 1 ,R(x)e~ iHtott . In order 
to compute $i ; /j(x, t) it is convenient to evaluate first 
bi(q,t) = e m *'* t bi{q)e~ iHtott in terms of the bi(q)'s and 
&j(g)'s, since they diagonalize Ho which defines the en- 
semble we average on. 

After a lengthy algebra involving a direct and an in- 
verse Bogolioubov transformation in Eq. and the in- 
verse of Eq. [21 one finds 



h(q,t) = A(q,t)h(q)+B*(q,t)bl(-q) 
+ C(q,t)b 2 (q)+D*(q,t)~bl(-q), 



(16) 
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where 

A(q,t) = coship(q) [— ismv s \q\tcosh2(p s (q) 
+ cosv s \q\t + (a «-> s)] 

— sinh</?(g) [—i sin v s \q\t smh2(p s (q) 
+(a <-> s)] , 

B*(q,t) = — sinhy>(g) [— i sin u s |g|t cosh 2ip s (q) 
+ cosv s \q\t + (a <-> s)] 
+ coship(q) [—i sxnv s \q\t sinh 2y> s (g) 
+(a <-> s)] , 
C(q,t) — cosh <p(q)[—i sin v s \q\tcosh.2<p s (q) 
+ cosv s \q\t — (a «-> s)] 

— sinh^g) [— isinu s |(7|isinh2(p ;s ((7) 
™(a <-> s)] , 

D*(q,t) = — sinh i^(g) [— i sin v s \q\t cosh 2ip s (q) 
+ cosv s \q\t — (a <-> s)] 
+ coship(q) [—isinv s \q\tsmh 2ip s (q) 
-(a «->*)]. (17) 

The above relations allow to evaluate exactly the time- 
dependence in the correlation function, which reads: 



. . 1 r / 2tt 
( x , t) = - — exp - 2^ —f 



x[(\A(q,t)\ 2 + \C(q,t)\ 2 )(-ismqx + 1 - cosga;) 
+ (\B(q,t)\ 2 + \D(q,t)\ 2 )(isinqx + 1 - cos gx)]](18) 

The sum over q is performed by using that q — 2nn/L 
and we end up with 

qW( x t ) = - 

R V ' ; M x \^(2cosh2p+cosb.(2tp-4<p !) )+cosh.(2<p—A(p a ,)) 



d(x — 2v s t)d(x + 2v s t) 



d(2v s t) 2 
d(x — 2v a t)d(x + 2v a t) 



^ (cosh(2<^> — 4<p B ) — cosh 2^p) 



d(2v a ty 



^ (cosh(2cp— 4ip a ) — cosh 2tp) 



(19) 



where d(x) = Lsin(nx/ L). In order to obtain the above 
expression we have introduced momentum-independent 
ip and ip Sta such that sinh /(g) « e - ^ ' 9 '/ 2 sinh /(0) with 
Rq of the order of the range of the e-e interactions; in 
this way the finite overall prefactor c is produced^}. 

Let us now consider some relevant limits. 

First, it is straightforward to check that in the non- 
interacting case (i.e. ip = ip s = <p a = and v a = v s — 
vp) we recover 



on the other hand for t 



G«(x,0) 



d(x) 
we find 



d(x) cosh2 v ' 



(20) 



(21) 



where d is again a finite constant; this result corresponds 
to the correlation function of an individual Luttinger liq- 
uid governed by H\ or H 2 , as it should. 

More interesting is the case in which g^' — g^' — 
g^ = (i.e. (p — and — ip a = ip s = Cp ^ and 
v a = v s = v Dp). In this case the system corresponds 
to have we two identical replicas of the model considered 
by Cazalilla, i.e. an isolated noninteracting ID system 
with density-density interactions switched on at t = 0. 
In this case the we obtain 

.(i) 



G^(x,t) 



rf(2;)3(cosh4 V +l) 



d(x - 2vt)d(x + 2vt) 



d(2vt) 2 



i(cosh4ip-l) 



(22) 



which recovers the result of Cazalilla (Eq. 9 of Ref. [l2|). 

Anyway the most remarkable result is obtained from 
the long-time limit t — > oo of Eq. In this case we get 



i (2 cosh 2ifi+cosh(2ip—4ip s )+cosh(2ip—4ip a )) 



X 4 



(23) 



where we have also taken L — > oo. We point out that 
if the infinite length limit is not taken, the correlations 
obey an oscillatory behavior even at long time, confirm- 
ing that in a finite-sized system a steady state is not 
reachable. However, as already discussed in Ref. 0, in 
mesoscopic systems with finite L (e.g. L ~ 1/xm for car- 
bon nanotubes) a genuine stationary state can be reached 
because of finite temperature effects. Indeed if the tem- 
perature T is larger than the level spacing Ae ~ vp/L, 
we have to replace L by vp/T and sin by sinh in the 
definition of d(x). As a consequence for very large dis- 
tances/times the Green function decays exponentially as 
G X R ~ e~' a: ' I/T , where v is the same exponent displayed 
in the power-law at T = 16] . Thus it appears that 
in the exponential regime we can still keep track of the 
Luttinger liquid correlations encoded by v. On the other 
hand for distances/times smaller than the thermal length 
v' F /T (with v'p the renormalized Fermi velocitv|l6j) the 
power-law discussed above still holds. 

In the thermodynamic limit, it is worth to notice that 
the result in Eq. 1231 does not coincide with the equi- 
librium equal-time correlation function of a two coupled 
Luttinger liquids governed by Hq + H\i. Indeed in the 
latter case the correlation function reads 



G^(x) 



| a ,|i(cosh(2y> a )+cosh(2^ Il )) ' 



(24) 



with finite c". The difference in the critical exponents 
in the two cases reveals that in the long-time limit the 
ground state of H evolved by H + 6(t)H\2 reaches a 
stationary state which is not the ground state of Hq + 
Hi2- This is due to the fact that for t > the energy is 
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conserved and therefore the system cannot relax to the 
ground state. Indeed the critical exponent in Eq. 1231 is 
larger than the critical exponent in Eq. [2U consistently 
with the fact that in the latter case the system is able to 
optimize the repulsive interactions given by H\2- 

Now let us discuss physical consequences of this result. 
To this end we introduce the critical exponent a govern- 
ing some relevant observables like the tunneling density 
of states o[ri| 



(25) 



which is detectable with scanning tunneling microscopy 
by probing the bulk the ID conductor[18] 19] . 

In order to compare to real systems, we have to intro- 
duce spin in the previous analysis. This is completely 
strightforward as long as we consider spin-independent 
e-e interactions. For spinless electrons, the a exponent 
is obtained by subtracting 1 to the exponent one-particle 
correlation function, while for the corresponding spinful 
system we just have to half such a value [2 5|. 

Anyway it is worth to note that, while in equilibrium 
systems the critical exponent of p(u>) can be safely ex- 
tracted from the equal-time one-particle Green function, 
in out-of-equilibrium situations such a procedure is far 
from obvious [20j. Therefore it would be convenient to 
relate a to some other energy-independent observable, 
like the momentum distribution, which also displays a 
power law|l7j]: 



n(q) <x sign(g - k F )\q - k F \ c 



(26) 



n(q) is accessible, for instance, from angle-resolved- 
photo-emission-spectroscpy, which has been success- 
f ully app lied to quasi-lD meterials[2"H [22j| to measure 
a 23} 24j. We also mention that such a technique is not 
affected by the uncertanties regarding the contacts, which 
are present in transport experiments. At this point, we 
note again that the above power-law would be cutoff by 
finite temperature effects if \q — kp\ ~ T jvp. 

From Eqs. 1231 and l2~4l the exponent governing n(q) in 
the two different cases discussed above read 



^asympt 
1 



(2 cosh 2ip + cosh(2(y9 — 4</? s ) + cosh(2</? — 4</? a )) 
o 

a quii = \(cosh(2ip s ) + cosh(2<y3 a )) - ^ . 



1 

" 2 ' 
(27) 



The model we have considered could find an experi- 
mental realization in ID fermionic systems built within 
optical potential. Unfortunately in these systems the in- 
teraction is short-ranged and this makes the term H12 
hard to realize. Anyway long-range dipole-dipole inter- 
actions have been recently obtained in chromium bosonic 
atoms [26) and Fermi gases with long-range interactions 
are likely becoming available. 



On the other hand it is tempting to relate our results 
to the physics of carbon nanotubes. Metallic carbon nan- 
otubes are believed to be rather good (although approx- 
imate) realizations of Luttinger liquids since in normal 
conditions the main correlation effects come from the 
long-range part of the Coulomb repulsion (through the 
g+ parametres), while the back-scattering interactions 
with large momentum transfer suffer a \l R sup pression, 
where R is the radius of the nanotube[18l||l9l|. Indeed 
back-scattering interactions are marginal operators in the 
renormalization group sense and in carbon nanotubes are 
predicted to breakdown the Luttinger liquid state only at 
very low energy scales p^j. This is supported by direct 
observation of power-law decay of the linear conductance 
as a function of temperature 27] 28] 29] . 

It is worth to recall that each individual carbon 
nanotube is composed itself by two coupled (identi- 
cal) Luttinger liquids, since there are a Left and a 
Right branch respectively around each of the two Fermi 
pointsHliill^- 

This implies that the critical expo- 
nent governing the decay of the correlation functions 
in these systems is actually one half the exponents de- 
fined in Eqs. [57| Therefore in the following we will use 



asympt 



aasympt/2 and a 



(nano) 
equil 



a e 



ail/2- 



The idea is to compare two possible measurements re- 
alized in two different experimental setups. The first sys- 
tem (setup A), which gives access to o4sym°pt, is made by 
two (initially distant) nanotubes which can be quickly 
aligned at controlled distance. This makes possible to 
mimic the switching on of inter-tube interactions. In 
this case, of course, the critical exponent could be only 
be measured in the final steady state and any transient 
regime is definitively unaccessible. The second system 
(setup B), which gives access to 0^™°^, consists in hav- 
ing the carbon nanotubes synthetized at a given distance, 
for instance in a double-wall nanotube. In this case the 
system has fixed inter-tube interactions built-in from the 
very beginning. 

Typical metallic nanotubes correspond to armchair 
(10,10) geometry (i.e. radius R ~ 7 A) with length 
L ~ 1/im. The strength of the density-density Coulomb 
interactions is then|l8j |30| 



<?i 2) (<7~o)= 5 n g ~o) 



^ 2) (g~o) = ^ 4) ( g ~o) 



4e 2 

KITVf 

4e 2 



-K (q c R) , 



K (q c D) , 



(28) 



where Kq is the modified Bessel function, q c = 2tt/L 
is the infrared momentum cutoff due to the finite size 
of the nanotubes, D > R is the intertube distance, 
k ~ 2 is the dielectric constant of typical nanotubes 
and vp as 8 x 10 5 m/s. We recall at this point that 
the amplitude of intra-tube backscattering interactions 
is approximativelv|l8j 0.1e 2 a/Kirv F R (with a = 2.46A) 
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FIG. 



1: Plot of the critical exponents a^sympt (dashed curve) 



(nano) 
equil 



(solid curve) for typical (10, 10) nanotubes as a 



and a 

function of the logarithm of the distance D between the nan- 
otubes. In log 10 D, D is expressed in A. 



which is much smaller of both intra- and inter-tube for- 
ward scattering interactions for D > R. We also mention 
that single-electron tunneling between nanotubes is usu- 
ally strongly suppressed since the misalignment of the 
respective carbon lattices prevents momentum conserva- 
tion in the tunneling process. Therfore the experimental 
setups that we propose could be fairly approximated by 
the model of two coupled Luttinger liquids. 

By inserting obtained in Eqs. in Eqs. and 
1101 we can estimate the critical exponents accessible to 
experiments within setups A and B. 

The two exponents are plotted in Fig. ^as a function 
of the distance between the nanotubes. For D — > oo the 



intertube couplings g^ and g^ ] vanish (see Eq. I28J) and 



(2) 



i (nano) (nano) ^ enc j ^ Q common value 



asympt 



equil 



(nano) 

v 

disc 



1 



(cosh 2^- 1) « 0.43. 



(29) 



which governs the correlations of two completely inde- 
pendent nanotubes, according to Eq. 1211 This value is in 
agreement with previous estimates given for single typi- 
cal nanotubes [19j. 

When the two nanotubes are brought close to each 
other in setup A (Fig. 1 dashed line), the intertube in- 



teraction are switched on and 



ticular at D 



R we find 

(nano) 
asympt 



^asympt 



0.60. 



increases. In par- 



(30) 



Regarding setup B (Fig. 1 solid line), if the two carbon 
nanotubes have been synthctized at the same distance 
D ~ R, the predicted exponent is 



(nano) 
oquil 



that is significantly smaller than ct^y^pt ■ We believe that 
such a large difference could be detected, even though the 
appropriate experimental setup might be difficult to re- 
alize. In particular in setup A the two nanotubes should 
approach each other at a speed comparable to (or even 
higher than) the Fermi velocity (~ 10 6 m/s). We note, by 
the way, that this problem does not arise in cold atomic 
systems, where the Fermi velocity is orders of magnitude 
smaller. 

In conclusion we have computed the time-dependent 
evolution of the single-particle correlation function of two 
Luttinger liquids coupled by a sudden switch-on of the in- 
ter liquid interaction. This allows to evaluate the critical 
exponent a governing some physical observables accessi- 
ble to real experiments, such as the momentum distribu- 
tion function. We find that the in the long-time limit the 
initial ground state relaxes to a stationary state which is 
not the equilibrium ground state of the coupled system. 
In the latter case the critical exponent a G q U ii results to be 
smaller than the asymptotic exponent a asym pt of about a 
factor 2. An experiment capable to detect such a remark- 
able finding in carbon nanotubes is proposed. Finally we 
believe that the present study could be also relevant in 
experiments involving ultracold fermionic atoms loaded 
in optical lattices, where tunable Luttinger liquids might 
be realized in the near future. 
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